Stochastic Semi-Classical Description of Fusion at Near-Barrier Energies 
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Fusion reactions of heavy ions are investigated by employing a simple stochastic semi-classical 
model which includes coupling between the relative motion and low frequency collective surface 
modes of colliding ions similarly to the quantal coupled-channels description. The quantal effect 
enters into the calculation through the initial zero-point fluctuations of the surface vibrations. A 
good agreement with results of coupled-channels calculations as well as experimental data is obtained 
for fusion cross sections of nickel isotopes. The internal excitations in non-fusing events as well as 
the fusion time are investigated. 
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I. INTRODUCTION 

Over the past years various models with different lev- 
els of approximation have been developed to explain the 
enhancement of sub-barrier fusion cross sections close to 
the Coulomb barrier energy. The description of fusion 
reactions is generally realized by phenomenological mod- 
els such as barrier penetration models (see ref. [l| for a 
review) or more sophisticated approaches like coupled- 
channels approach These models provide some 
qualitative as well as quantitative explanations to the 
enhancement of the cross sections as a result of coupling 
to selective collective modes. In such macroscopic ap- 
proaches, one of the main ingredients is the nuclear inter- 
action potential which can be taken as the Bass potential 
@ , the proximity potential @, 0] , or the double- folding 
potential Q. Some microscopic approaches equipped 
with full quantum treatment such as Time Dependent 
Hartree-Fock model are also used to give quantitative ex- 
planations of nuclear fusion reactions by providing a con- 
nection between the macroscopic and microscopic phe- 
nomena 

In this paper, we investigate fusion reactions of heavy- 
ions by employing a simple stochastic semi-classical 
model. In this model, the coupling between the rela- 
tive motion and low frequency collective surface modes 
is incorporated into the description in a manner similar 
to quantal coupled-channels calculations. Quantal effects 
enter into the calculation through the initial zero-point 
fluctuations of the surface modes. In the calculations, 
we include the coupling of relative motion with the low- 
lying quadrupole and octupole surface modes. The zero- 
point fluctuations of surface vibrations are simulated in 
a stochastic approximation by generating an ensemble 
of trajectories in accordance with the Gaussian distribu- 
tion of zero-point fluctuations. The zero-point quantum 
fluctuations of surface modes lead to barrier fluctuations, 
which enhance the fusion cross-section at near and sub- 
barrier energies. This stochastic semi-classical model has 
been proposed in ref. [l4| (see also [l5| for a detailed re- 



view). However, only a few applications of the model 
have been carried out so far |16l4l8l |. This is probably 
due to insufficient computation power at the time it was 
proposed. In this work, employing this stochastic model, 
we investigate fusion reactions of nickel isotopes at near- 
barrier and sub-barrier energies. We compare our re- 
sults with the experimental data as well as the quantal 
coupled-channels calculations with the same input pa- 
rameters. In section two, the model is briefly described. 
In section three, calculations of the fusion cross-sections 
of nickel isotopes are presented. In section four, further 
applications of the approach are illustrated and conclu- 
sions are given in section five. 



II. STOCHASTIC SEMI-CLASSICAL MODEL 

In order to describe heavy-ion collisions and fusion at 
near-barrier and sub-barrier energies, we follow the idea 
originally proposed by Esbensen et al. [l4| • In this semi- 
classical model, heavy-ion collisions are described by a 
Hamiltonian in which the relative motion is coupled to a 
number of low-lying collective modes (surface vibrations) 
of the colliding ions. Treating the surface vibrations in 
harmonic approximation, the semi-classical Hamiltonian 
for the colliding ions is given as, 
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where R represents the relative distance between two cen- 
ters of mass of the colliding nuclei and P is the corre- 
sponding relative momentum. In this expression, the first 
and second terms are the radial kinetic energy and the 
rotational kinetic energy with orbital angular momen- 
tum I. The quantities Vc(R) and Vn(R, repre- 
sent the Coulomb potential energy and the nuclear inter- 
action potential. The parameters set £1 = {^1,02,^3} 
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in the nuclear potential describes rotation angles of the 
vibration axes of the nuclei, which are specified below. 
The last term in Eq. (TTJ) is the Hamiltonian for 2N har- 
monic oscillators corresponding to the vibrational modes 
(A = 0, N — 1) of projectile and target ions (i = 1, 2). 
The quantities au, IT a, -DiA, and Cu indicate variables 
for the vibrational modes, that correspond to the defor- 
mation variables, the corresponding momenta, inertia pa- 
rameters, and spring constants, respectively. The spring 
constants and the inertia parameters of harmonic oscilla- 
tors are determined in terms of the deformation param- 
eters and the excitation energies E* x of the modes 
according toCa= and D iX = h 2 /2E* x (3 2 x , re- 

spectively [l8|, LL9| ■ In the ground state, the variances of 
vibrational variables and the variances of the correspond- 
ing momentum are expressed according to cr a . x = (3 iX and 

°"n ix = h/20i\- 

The classical equations of motion for the relative dis- 
tance and the vibrational variables are given by, 
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In the model, a dissipation of the relative energy occurs 
due to excitations of the surface modes. An important 
dissipation mechanism due to nucleon exchange between 
projectile and target nuclei 0, HH is neglected here. 
However, at low bombarding energies, the mechanisms 
related to excitations of the surface modes dominate com- 
pared to the nucleon exchange mechanism (see Figure 1 
of rcf. [22]). Consequently, the model provides a deter- 
ministic description of the average properties of collision 
dynamics at low energies. 

Due to a short de Broglie wavelength, the classical ap- 
proximation works well for relative motion at near-barrier 
energies where the effect of tunneling is small compared 
to that of the surface excitations. However, as a result 
of a few phonon excitations during the collision process, 
the dynamics of surface vibrations is far from the clas- 
sical limit, and should be treated in a quantal frame- 
work. A standard description is provided by the quantal 
coupled-channels calculations. Here, instead of the stan- 
dard coupled-channels description, we include the quan- 
tal fluctuations of the surface modes by incorporating the 
initial zero-point fluctuations of the vibrational modes. 
We can determine the phase space distribution function 
F(a, n) of a harmonic oscillator ground state, by taking 
the Wigner transform of the ground state wave function 
to find, 
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where a = ai\ and n = Tli\ are the variances of coor- 
dinate and momentum distributions for each vibrational 
mode. The quantal zero-point fluctuations of the vibra- 
tional modes are incorporated in a stochastic manner. An 
ensemble of trajectories is generated by solving the clas- 
sical equations of motion with initial conditions ai\(0) 
and n^O), that are randomly selected according to the 
corresponding Wigner distribution F(a,U). Once the en- 
semble of trajectories is generated, different observables 
are calculated by averaging over the ensemble. 

The following approximation for the Coulomb poten- 
tial is employed, 
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where Rc = Ri + R% is the sum of the equivalent sharp 
radii [23, 24]. The nuclear part of the interaction is com- 



puted using the double-folding potential as, 



xV NN (R -fi+ r 2 )d 3 r 1 d 3 r 2 . (5) 

Fig. 1 illustrates the geometry of two colliding ions in 
which xi,yi,zi and x 2 ,y 2 ,z 2 denote two sets of coor- 
dinate systems with fixed orientations and origins that 
are attached to the centers of the nuclei. The red lines 
indicate the vibration directions. The angles between po- 
sition vectors f*i, r 2 and the axes z\^ 2 are indicated by 8\ 
and #2- The angles Q' x and 9' 2 represent the angles be- 
tween position vectors f\ , r 2 and the vibration directions 
of the nuclei, respectively. The nuclear folding potential, 
in addition to the relative position R and the vibrational 
variables ctix, also depends on three independent angles 
£1 = {f2i, ^2, ^3} which specify the relative orientation of 
vibration directions of the nuclei. The angles f2i and Q 2 




FIG. 1: (color online) A schematic view of nuclei with 
quadrupole vibrations. The red lines indicate the vibration 
directions. 

are defined as the rotation angles of the vibration direc- 
tions about xi and x 2 axes, respectively. These two an- 
gles define orientations within the same plane {x\. 2 = 0), 
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and in order to cover all possible orientation configura- 
tions, one more angle f2 3 , is needed to rotate the oscilla- 
tion axis of one of the nuclei about the corresponding y or 
z axis to account for the off-plane orientations. For con- 
venience, a rotation about the y axis is considered. The 
rotation is not indicated in the figure. We approxi- 
mate the nuclear density distributions by two-parameter 
Fermi functions, which can be conveniently expressed in 
terms of the angles 8[ and 8' 2 as 



POi 



1 + exp [(n - Ri{8' i ,aix))/a i 



where poi is a normalization constant, is the diffuseness 
parameter and Ri(8' il aix) denotes the deformed nuclear 
radius of each nucleus. For small amplitude vibrations, 
this quantity is expanded in terms of spherical harmonics, 

= R 0l ^l + ^aa^^PA(co S ^)j(7) 

From the geometry displayed in Fig. 1, it is possible to 
deduce the following relations between angles, 



cosf^ = cos fli cos 9i — sin fix sin 8\ sin <j>i, 



(8) 



cos 8' 2 — cos Vt 2 (cos cos 62 — sin O3 sin 82 cos 02 ) 

— sin 5^2 sin ^2 sin 02 ■ (9) 

The normalization constants poi in Eq. (|6]) are obtained 
from the equation, 



pi(ri,Qi,aix)d 3 ri = A4 



(10) 



where A, is the mass number of the i th nucleus. For the 
equivalent sharp radii of the spherical nuclei, we take the 

1 /3 

values given by R 0l = 1.31AV -0.84 fm. 

A global description of the nucleus-nucleus potential 
can be achieved via the folding potential Eq. (JS|) by con- 
sidering a zero-range nucleon-nucleon interaction, 



V NN (r) = V S(r), 



(11) 



which is equivalent to a finite-range nucleon-nucleon in- 
teraction |25l - l28j |. Then, the nuclear potential Eq. ([5]) 
becomes 

V N (R,Q,aix) = 

V / p 1 (r 1 ,9' 1 ,a 1 x)p2(r2,O2,a2x)rldrid(cos0 1 )d(l) 1 ,(12) 



where 8[ and 9' 2 are given by Eq. ([8]) and Eq. ([9]), re- 
spectively with 



cos 82 = (ri cos 81 — R)/r 2 , 
r2 = 



(13) 



\ + R 2 -2r 1 Rcosd 1 , (14) 

(15) 



In order to simplify the numerical simulations, a further 
approximation is introduced. For small amplitude vibra- 
tions, Taylor expansions of the nuclear densities and the 
nuclear potential are introduced to the first order around 
(*ix = to give, 
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respectively. 

Fig. 2 shows examples of potential energies in head-on 
collisions of two 64 Ni nuclei for different orientations and 
deformations as a function of the relative distance. In 
this figure, for simplicity, only the quadrupole vibrations 
of both nuclei in two different orientations are shown. 
One of them corresponds to the case when vibrations are 
along the direction of the relative motion indicated by 
blue curves. In the following, this orientation is referred 




FIG. 2: (color online) The s-wave potential barriers of 64 Ni 
+ 64 Ni system are plotted versus center to center distance 
for different quadrupole deformations and orientations. The 
blue curves represent the oblate (upper curve) and prolate 
(lower curve) deformations along the relative motion. The 
red curves stand for the prolate (upper curve) and oblate 
(lower curve) deformations that are perpendicular to the rel- 
ative motion direction. The bare potential with no defor- 
mation is indicated by the black color. The potentials are 
computed by summing the Coulomb Vc and nuclear Vjv con- 
tributions with deformation variables equal to their variances, 
jQi2| = \ct22 \ = P2 = 0.215. The zero-range potential strength 
in Eq. (TlT]) is taken as V = -456 MeVfm 3 



to as the ZZ configuration. The other orientation cor- 
responds to the case when vibrations are perpendicular 
to the direction of the relative motion indicated by red 
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curves. We refer to this orientation as the YY configura- 
tion. It is observed that in the ZZ configuration, the dif- 
ference between the potential barriers for oblate and pro- 
late deformations is larger than the difference in the YY 
configuration. Any orientation in between these two con- 
figurations leads to a difference in barrier heights, which 
is smaller than that of the ZZ configuration and larger 
than that of the YY configuration. Hence, the ZZ and 
YY configurations represent the extreme states of bar- 
rier fluctuations due to surface vibrations. When the off- 
plane orientations are included (Q3 ^ 0), the difference in 
barrier fluctuations becomes minimum for the configura- 
tion where vibrations are perpendicular to the direction 
of the relative motion and perpendicular to each other. 
This orientation is referred to as the XY configuration, 
which is not included in Fig. 2. It is clear that different 
orientations of the nuclear surface vibrations have a very 
large influence on the fusion barrier fluctuations. There- 
fore, we need to calculate any observables by averaging 
over all possible relative orientations of the vibration di- 
rections. We carry out this averaging by sampling all 
three angles, {f2i, ^2, ^3}, from a uniform distribution 
in the interval [0,27r]. 

III. FUSION CROSS-SECTIONS 

Coupled-channels calculations are often employed for 
describing fusion cross-sections at sub-barrier energies. 
These investigations indicate that low-lying surface 
modes such as 2 + and 3~ make the dominant contri- 
bution to sub-barrier cross-sections [29l - l3Tj |. Retaining 
only these two modes, we carry out stochastic simula- 
tions to describe the fusion process of Nickel isotopes. 
In order to compare our results with that of coupled- 
channels calculations of Nobre et al. ;32|' , we adopt the 
same parameters as in that reference. It is important to 
note that none of these parameters are adjustable. The 
quadrupole (A = 2) and octupole (A = 3) deformation 
parameters are fa = 0.215, fa — 0.263 for 64 Ni, and 
fa = 0.205, fa = 0.235 for 58 Ni. The excitation ener- 
gies are E% = 1.35 MeV, E£ = 3.56 MeV for 64 Ni, and 
E% = 1.45 MeV, E% = 4.48 MeV for 58 Ni. The zero- 
range potential strength and the diffuseness parameter 
are Vq = —456 MeVfm 3 and a = 0.56, respectively [25| . 
For simplicity, we consider that the quadrupole and oc- 
tupole vibrations of each nucleus are aligned in the same 
direction. 

We calculate the fusion cross section using the stan- 
dard expression, 

t2 Imaas 

^ 1=0 

where E and Pi(E) represent the incident center of 
mass bombarding energy and the partial transmission 
probabilities, respectively. Fig. 3 shows the results 
of the stochastic calculations and comparisons with the 
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FIG. 3: Fusion cross sections of 64 ' 58 Ni + 64 ' 58 Ni systems cal- 
culated with the coupled-channels model (CC) and the Frozen 
Approximation model (FA) of Nobre et al. [321 ] are compared 
with our results (the solid line). Barrier heights are indicated 
in the figures. The shaded areas are explained in the text. 



coupled-channels calculations. The solid line corresponds 
to the cross section obtained by averaging over all orien- 
tations. The shaded area in the figure illustrates the 
cross section fluctuations due to the effect of different 
vibration orientations. The upper (lower) boundary of 
the shaded area corresponds to the cross section when 
vibrations are along (perpendicular to) the relative mo- 
tion direction which gives rise to maximum (minimum) 
cross-section due to the largest (smallest) barrier fluc- 
tuations. The coupling to surface vibrations increases 
the fusion cross-sections at sub-barrier energies, whereas 
it decreases the cross-sections at over-barrier energies. 
This effect is further enhanced due to the orientation 
configurations which can be easily visualized by looking 
at Fig. 3. At very low bombarding energies, the stochas- 
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tic fusion cross-sections tend to approach the ZZ con- 
figuration limit since the barrier fluctuations are larger 
in this case. On the other hand, at over-barrier ener- 
gies, the fusion cross-sections tend to approach the XY 
configuration limit since the transmission is reduced by 
very large barrier heights appearing in the ZZ configu- 
ration. The reduction in the ZZ configuration is eventu- 
ally increasing the weight of events in the XY configura- 
tion. In this figure, we also compared our results with 
the coupled-channels calculations (CC) and the frozen 
approximation (FA) of Nobre et al. !32] which is an ap- 
proximation to the coupled-channels calculations. Even 
though same parameters are used in both calculations, 
there are some minor differences. While they used the 
double folding potential in their CC calculations, a frozen 
density approximation is further assumed. Both calcu- 
lations, CC and FA, are performed with a parabolic ap- 
proximation for the potential barriers with effective cur- 
vatures. The stochastic approach does not have these 
drawbacks, which are important at low sub-barrier en- 
ergies. Furthermore, in the stochastic description, the 
effects of vibration direction orientation of the nuclei is 
incorporated in a natural way. On the other hand, the 
stochastic model does not include the effect of quantum 
tunneling. Consequently, the stochastic approach pro- 
vides a good description at near-barrier energies where 
barrier fluctuations due to the surface vibrations provide 
dominant contributions to the transmission probability. 
In Fig. 3, it is seen that the stochastic semi-classical 
description provides a good approximation to the quan- 
tal CC model and its approximate FA version for nickel 
isotopes fusion at near-barrier energies. In Fig. 4, the fu- 
sion excitation functions obtained by the stochastic semi- 
classical model are compared with three different data 
sets. Our results are in good agreement with the data 
near and below the barrier. The model over-estimates 
the experimental cross-sections at over-barrier energies 
due to the following reasons. First, the nuclear potential 
that we employ is energy independent. Energy depen- 
dence of nuclear potential can have a significant effect 
on fusion mechanism above the barrier. Actually, there 
are contradictory claims for the energy dependence of 
nuclear potentials. In ref. [25l427j . it is claimed that, 
at near-barrier energies, the energy dependence vanishes 
while at higher energies it becomes important. On the 
other hand, in some microscopic calculations, an oppo- 
site behavior is found [l^. Second, the folding potential 
that is employed here has a deep minimum as a result 
of the zero-range interaction. It has been shown that a 
repulsive core within the M3Y nucleon-nucleon interac- 
tion give s rise to shallower, more realistic nuclear poten- 
tials [36[ . The maximum angular momentum l max , up to 
which the partial transmissions are summed over in Eq. 
(|18p. is determined by the fact that the potential pocket 
disappears. For deep potentials, l max is larger than that 
for shallow potentials, which increases the over-barrier 
cross sections. We also believe that the lack of dissipa- 
tion in the relative motion as well as in surface vibrations 



is the third reason for the over-estimation of data above 
the barrier. Dissipation mechanism can be easily incor- 
porated into the stochastic description. Nevertheless, it 
is not considered in this study. Our task in this work is to 
test the stochastic semi-classical model by carrying out 
simulations of fusion mechanism and compare the results 
with the quantal coupled-channels calculations. There- 
fore, we avoid employing a complicated nuclear potential. 
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FIG. 4: The fusion cross sections of 64 ' 58 Ni + 64 ' 58 Ni systems 
calculated with the stochastic zero-point model are compared 
with the data of ref. .'!•"> . Barrier heights are indicated in 
the figures. 



IV. DISCUSSIONS 

The stochastic semi-classical model provides a simple 
framework not only for fusion cross-sections but also to 
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evaluate some relevant observables such as the time dis- 
tribution for fusion and the kinetic energy distribution 
of non-fusion events which can hardly be accessed in a 
fully quantal framework. Here, we restrict the study to 
collisions between 64 Ni nuclei. Fig. 5 shows five illus- 
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FIG. 5: Fhe center of mass energies for five sample events, 
with incident energy of 90 MeV, are plotted versus the relative 
distance of two 64 Ni nuclei. Fhe center of mass energy is the 
sum of the first four terms in Eq. {TJ, hence it is the energy 
of the relative motion. Fhe angular momentum is set to zero. 
Fhe s-wave potential barrier is shown for the bare case where 
there is no coupling to the surface modes (thick line). 

trative events as a function of relative distance of 64 Ni 
ions for head-on collisions at bombarding energy of 90 
MeV. In this figure, the thick line indicates the bare po- 
tential. Even though the bombarding energy is below 
the bare barrier, as a result of barrier fluctuations, three 
of the five events end up fusing, while two events after 
inelastic collision re-separate. During inelastic collisions, 
part of the incident energy is dissipated by excitations 
of the surface vibrations. Using the ensemble of events 
generated for description of the collision process, we can 
calculate the final kinetic energy distributions of the non- 
fusion events. The upper panel (a) of Fig. 6 shows final 
kinetic energy distributions of non-fusion events at dif- 
ferent incident energies in head-on collisions of 64 Ni + 
64 Ni systems, while the lower panel (b) presents the final 
kinetic energy distributions for different relative angu- 
lar momenta at an incident energy of 96 MeV. As it is 
expected, the mean final energy increases with the inci- 
dent energy as well as with the angular momentum. Fhe 
variances of the distributions increase as incident energy 
increases and angular momentum decreases. Indeed, in 
this case, the rates of non-fusion events decrease eventu- 
ally increasing the uncertainty of the final energy. Fig. 
7 illustrates the average value of the final energy of the 
relative motion of non-fusion events. It is observed that, 
at very low incident bombarding energies, almost no en- 
ergy is transferred to surface excitations. As the energy 
increases, the amount of dissipated energy into the sur- 
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FIG. 6: Fhe distribution of the final energy of non- fusing 
events in collisions of 64 Ni nuclei is plotted for different in- 
cident center of mass energies (a) and different angular mo- 
menta (b). For each plot the distribution function is normal- 
ized with the total number of scattered (no-fusion) events. 



face modes increases. After the energy exceeds the fusion 
barrier height at Vb — 96.5 MeV, the mean final energy 
is almost constant indicating that the total excitation 
energy E* is almost linearly increasing with the incident 
energy. In the simulations of the trajectories, we take 
the initial relative distance as R — 20 fm. We call the 
event as a fusion event, if the trajectory evolves all the 
way until the separation distance reaches to R = 5 fm. 
For fusion events, we define the time it takes to travel 
from R — 20 fm to R = 5 fm as the fusion time. Fhe 
upper panel (a) of Fig. 8 indicates the distribution of 
fusion times at different incident energies in central col- 
lisions of nickel ions, while the lower panel (b) shows the 
fusion time distributions for three different orbital angu- 
lar momenta at incident energy of 96 MeV. Again, as we 
expect, the mean fusion time increases with increasing 
angular momentum and decreasing incident energy. 



V. CONCLUSION 

Employing a stochastic semi-classical model for low- 
energy heavy-ion collisions, which was proposed origi- 
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FIG. 7: The average values of the total excitation energy (a) 
and the final energy (b) in head-on collisions of 64 Ni nuclei 
are plotted as a function of the incident center of mass energy. 
The excitation energy is given by E* = Ei nc — Ea na ,i- 



nally by Esbensen et al. [14] . we carry out simula- 
tions for describing fusion and some other properties of 
the collision of nickel isotopes at sub-barrier and near- 
barrier energies. As known from coupled-channels calcu- 
lations, dominant effects for describing sub-barrier fusion 
arise from the coupling of the relative motion with the 
low-lying surface vibrations. Therefore, in analogy to 
coupled-channels calculations, we consider the coupling 
of the relative motion with low-lying collective surface vi- 
brations. Since the de Broglie wavelength is very short, 
a classical treatment provides a good approximation for 
the relative motion. However, the qauntal aspects of sur- 
face vibrations play a dominant role in the sub-barrier 
fusion mechanism. In the present approach, we incor- 
porate quantal zero-point fluctuations of the surface vi- 
brations in a stochastic approximation. In the applica- 
tions presented here, only the quadrupole and octupole 
vibrations of the projectile and target nickel ions are in- 
cluded. An ensemble of trajectories of the relative motion 
are generated by picking the initial conditions of surface 
vibrations according to the quantal zero-point fluctua- 
tions of their ground states. Some observables such as 
fusion cross-sections and final kinetic energy distributions 



of non-fusion inelastic collisions have been estimated by 
averaging over the ensemble of trajectories. The simple 
stochastic semi-classical approach provides a surprisingly 
good agreement with quantal coupled-channels calcula- 
tions as well as the experimental data for fusion cross- 
sections of nickel isotopes at sub-barrier and near-barrier 
energies. Of course, barrier penetration is not included 
in the description. However, at near-barrier energies, the 
dominant effects on fusion arise from the barrier fluctu- 
ations, which are well accounted for by the stochastic 
approach. 
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FIG. 8: The distributions of the fusion time in collisions of 
64 Ni nuclei are plotted for different incident center of mass 
energies (a) and different angular momenta (b). For each plot 
the distribution function is normalized with the total number 
of fusion events. 
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